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Abstract 

A new equation for the density contrast is derived in the 
framework of reexamined Newtonian cosmology taken into ac- 
count adiabatic matter creation in the universe. 

The age of the universe and the reach of non linear regime 
of the density contrast are usually treated separately in the 
literature and this may lead to controversial conclusions re- 
garding the most adequate scenario to describe the universe. 

We relate the age of the universe and the growing mode of 
the density contrast by introducing a variable £ that relates 
both of them, thus both aspects are treated simultaneously. 
We apply this procedure to the Friedmann type model where 
the source of particle production is ^ = 3nj3H. 



1 Introduction 



In 1934, E.A.Milne analyzed the expansion of the universe using el- 
ementary Newtonian theory (Tp. In the same year W.H.Mc Crea and 
E.A.Milne demonstrated the identity between the governing differ- 
ential equations for relativistic and Newtonian cosmology [@] . Three 
decades after, E.R.Harrison (3| extended Newtonian cosmology to 
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include pressure generalizing Mac Crea and Milne's approach for 
dust matter. The Harrison's approach is based on the following set 
of equations: 



dp 
di 




(1) 



(2) 



which are respectively the continuity equation, Poisson's equation 
and the motion equation. Posteriorly, Harrison [§] derived the mo- 
tion equation by introducing a non uniform pressure, namely: 



In the above equations p, P, u and denote respectively the density 
of matter, the thermodynamic pressure, the velocity field and the 
gravitational potential. 

The field equations for a homogeneous and isotropic universe 
are obtained by using equations (1), (2) and (3). These equations 
are identicals to the ones obtained in the framework of general 
relativity in the comoving coordinate system. 

However, it was observed in the end of the appendix II of the 
R.K.Sachs and A.M Wolfe work that the perturbed modified field 
equations to the first order are not the same in the general rela- 
tivity and Newtonian cosmology. Therefore, the density contrast 
equation for both theories are quite different. 

This ambiguity was solved by Lima, Zanchin and Brandenberger 
|6p that derived the continuity equation in the modified Newtonian 
approach to cosmology when pressure effects are included. The 
new continuity equation is completely consistent with relativistic 
perturbation theory and solves the contradiction pointed out by 
Sachs and Wolfe. 

In this paper, we generalize the evolution equation for the den- 
sity contrast obtained by Lima et. al., for the universe with particle 
production and obtained the growing and decaying modes for the 
density contrast. 



(JUL ~~ * ~~ * 

— + U ■ V r U = -V r <j) - (p + P th ) VrPth ■ 



(3) 
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In our work we introduce a variable (£) that depends of the age 
of the universe and of the growing mode of the density contrast. 
The use of £ can discern the cosmological model that furnishes a 
better concordance with the experimental evidence about the age 
of the universe and the time of structure formation 

This paper is organized as follows: we obtain in section II, in 
the linear approximation, the time evolution equation for the mass 
density contrast for a Newtonian universe with matter creation. In 
section III we derive the time dependence of the density contrast 
in which the source of particle production is ^ = 3n(3H. We ob- 
tain a range for the parameter (3 that updates the range obtained 
in the literature when the question of the age age of the universe is 
involved. Finally we use the reach of the non linear regime to es- 
tablish a criterion about the most adequate cosmological scenario. 



2 Equation for the density contrast 

In the context of matter creation the thermodyamical conservation 
law reads: 

dp 1 , „ p \ h d , T T . 

where n is the particle number density, h = p + P th is the enthalpy 
per unity volume and P th is the thermodynamical pressure. Con- 
sequently the continuity equation becomes 

! + v r .(pig + ^v r .tr=£*, (5) 

where ^ is the source of particle production. The other hydro- 
dynamical equations that describe the cosmic fluid are the momen- 
tum conservation equation and Poisson's equation, respectively: 

8ii P 

+ ^ r u = -V r $ - (p + -VVrP , (6) 

at c z 

and 

IP 

V r 2 $ = 47rG(p + — ), (7) 
cr 
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$ is the gravitational potential and P is the total pressure that 
includes the thermodynamical pressure and the creation pressure 
P , so 

P = P + P th . (8) 

The creation process will be considered adiabatic, consequently the 
pressure creation is given by [|7| : 

In order to obtain a cosmological scenario with matter creation one 
still needs to provide the state equation 

P th = vp- (10) 

Using the field equations for FRW model with particle pro- 
duction and the state equation above, we can write a differential 
equation for the scale factor: 

v 2 2nH ' V 2 2nH J 

Our main aim in this section is to find the differential equation 
that governs the evolution of the density contrast in a Universe 
where the particle production is given by ^ = 3n(3H, where (3 is 
a constant. Note that j3 = implies ^/ = and we recover the 
standard Friedmann model. 

To study the evolution of small fluctuations in an expanding 
universe we consider the standard perturbation ansatz: 

p = Pb (t){l + 6(r,t)}, (12) 

P th = Pthb(t) + 5P th (r,t), (13) 

$ = $ b (r,t) + <p(r,t), (14) 

u = u b (t) + v(r,t), (15) 

P = P b (t) + 5P(r,t). (16) 

The quantities carrying the subscript b represent the homoge- 
neous solution to the unperturbed equations. Inserting the above 
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expressions into eqs.(2),(3) and (4) we get to first order in perturba- 
tions 

. R v-p(u + l) vs 

v + r = -R Vv - ( V+ i)(i-m -R 1171 

VV = ^Gp b (t)R 2 5{l + 3u- 3(i/ + l)f3} (18) 
5 = ^Vv(l-{3). (19) 

ri 

Change to comoving coordinates following standard lines [|§ , elimi- 
nating the peculiar velocity from equations (17) and (19) and using 
(18), is ready that 

6 +2^5 + AnGp b (t)6{l + 3u- 3{u + 1)(3}{v + l)(/3 - 1) 

= (u-P(v+l))^. (20) 

This is the differential equation that governs the evolution of den- 
sity contrast in the presence of matter creation, when we describe a 
matter distribution with uniform pressure using the modified New- 
tonian equations. It is spatially homogeneous, so we expect to find 
plane wave solutions. The difficulty to apply Newtonian cosmology 
in the study of scalar perturbations is related with the impossibil- 
ity of the Newtonian scenario to describe long wavelengths pertur- 
bations. The cosmological perturbations can be of the two kind 
according to the wavelength A, so that A > d H or A < d H , where dh 
is the Hubble sphere. Perturbations with wavelengths larger than 
Hubble scale requires some form of a general relativistic theory 
of perturbations. For wavelengths smaller than Hubble scale the 
evolution of mass density can be studied using Newtonian theory. 
Application of Newtonian equations is further restricted to non- 
relativistic matter and cannot be used for relativistic components 
even for scales smaller than Hubble radius. 

The reexamined Newtonian equations indicate a way of obtain- 
ing the same evolution equation for the density contrast as could 
be by the full relativistic approach. In this way one can extend the 
domain of validity of equation (20) in order to analyze perturba- 
tions even in the radiation dominated phase and to apply the large 
wavelength limit. 
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Next, we derive the solution of the density contrast when the 
source of particle production is \P = 3n/3H, concentrating our at- 
tention on the growing mode. The principal advantage of this ap- 
proach is the relative facility to obtain the perturbed equations for 
the density contrast which is the same equation as in general rela- 
tivity when (5 = in the synchronous and comoving coordinates. 

Although the gauge invariant approach is conceptually more at- 
tractive since there is no need to identify the physical and unre- 
alistic space-times, nevertheless it is more complicated and the 
physical meaning of gauge independent variables do not in general 
possess any simple interpretation. Besides, the gauge mode that 
emerges from the differential equation for the density contrast is 
generally related with some kind of decaying mode when the uni- 
verse is homogeneous and isotropic |TTJ . 

3 Choosing a model 

If the space time has uniform curvature, the line element is given 

ds 2 = dt 2 - R 2 (t)( dr i - + r 2 d6 2 + r 2 sin 2 (e)d(j) 2 ) . (21) 
1 — kr z 

The evolution of the scale factor is obtained from solutions of dif- 
ferential equation [0 

RR + A(k + R 2 ) = , (22) 

where 

A = ^(z/+l)(l-/3)-l. (23) 

Explicit solutions for differential equation (22) when k = ±1 for any 
value of p are not easy to obtain. We have managed to find a class 
of exact solutions for the following cases: 

K = (3 1 — > R = R (l)^ (24) 

to 

K = (3=1 — ► R = R oe i , (25) 
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K=±1 , = 1 _K=i { ^l±%. ,26, 

1 e ±c ?(t + C 2 )CZ 

The subscript alludes to the present time and C\ and C 2 are inte- 
gration constants. 

Whether the Universe turns out to be spatially closed, open or 
precisely flat remains an empirical question. The recent results of 
the Boomerang provide convincing evidence in favor of the stan- 
dard paradigm: The Universe is flat, however the possibility that 
the universe might be spatially closed do not be discarded 
Consequently, we consider the Universe as flat for subsequent cal- 
culation. 

Using the scale factor (24) then the mass-energy density is given 

by: 

feG "'=(ATW (27) 

Using equations (27), (24) and the long wavelength limit into equa- 
tion (20) we finally obtain 

S + j5+(a-2)^ = 0, (28) 



where 



° = 3(, + l)(l-fl - 129 ' 



Solutions for equation (28) are 

s_ = c l r 1 , (30) 

S + = C 2 t 2 ~^-w«+v . (31) 

In the absence of particle production, /3 = 0, relations (30) and (31) 
become the usual density contrast for Friedmann models without 
creation Although the decreasing mode can be important in 



some circumstances, the increasing mode is responsible for the 
formation of cosmic structures in the gravitational instability pic- 
ture. Taking into account the decreasing mode the universe would 
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not have been homogeneous in the past. Besides, Peebles argued 
that a growing mode that starts to grow just after the end of radia- 
tion era has a negligible component of decaying solution [[T5] . 

One of the problems of the standard FRW model is related to 
the age of the universe; it is younger than it should be according 
to the experimental values. However, in the framework of particle 
production the value of the age of the universe is more compati- 
ble with the experimental evidences. To show this we define the 
quotient ( given by: 

c=|r, (32) 

at the same time coordinate. The subscript F and C refer respec- 
tively to the Friedmann standard model and the model with particle 
creation. 

Considering the scale factor (24) one obtains 

C = T^- (331 

If we expect that the open system model furnishes an oldest uni- 
verse, then C, > 1 implying < /3 < 1. Using the oldest globular 
clusters, B. Chaboyer [[T3J infers a range for the age of the universe, 
namely 9.6Gyr < t universe < lh.AGyr. Consequently, the Hubble con- 
stant for the standard model lies in the interval 

43.3 < H < 69.4 . (34) 

This range does not agree with the interval for H estimated re- 
cently by Willick and Batra p~5| ] , namely: 



80 < H < 90 . (35) 

One can easily verify that there is a small overlap between (34) and 
the values obtained by Richtler and Drenkhahn [ |16H , given by 

68 < H < 76 . (36) 

Using the scale factor (24) for dust matter universe, the Hubble 
function for a universe with particle production is given by 
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We can conciliate the ranges (35) and (36) using the Hubble func- 
tion above, where the creation parameter (3 lies in the interval 

0.13 < < 0.52. (38) 

This result updates the interval for (3 determined by Lima et al (SJ 
(0.34 < (3 < 0.60). 

Now, we study the reach of non linear regime, identifying the 
epoch for the reach of the non linear regime with the epoch for the 
formation of the super clusters. In addition, we suppose that the 
star formation is posterior to the formation of the super clusters. 
First, because the density of the super cluster in the universe is 
approximately equal to the universe density. Second, considering 
the universe to form "from top down" is acceptable for some type 
of models, for example the hot dark matter model [ |17| ]. In other 
words, we want to infer that the star formation is posterior to the 
super cluster formation and, consequently, to the reach of the non 
linear regime. The growing mode for the density contrast for the 
FRW standard model can be written as 

8 + = 5 d {U 2 '\ (39) 
td 

where the subscript d refers to decoupling time. Typically we have 
)dec < (10 -2 - 10" 3 ) ||17J, for a decoupling time of the order 10 5 ys . 
Substituting the decoupling time and the corresponding anisotropy, 
the reach for the non linear regime for FRW standard model occurs 
in the range 10 _1 Gyr < t < 3.16 x Gyr. Note that, the standard 
model furnishes a sub estimate age for the universe, around 8 
Gyr. Now, if we consider the difference between the age of the uni- 
verse and the reach of the non linear regime we obtain the interval 
5Gyr <t< 7 Gyr for the age of the super cluster. In the literature it 
is known that the age for the oldest stars in our galaxy is around 
12 Gyr 

The open system cosmology furnishes a reasonable solutions for 
the age of the universe. What happens with the reach of the non 
linear regime for the density contrast? 

To solve this questions we introduce a variable £ that relates the 
age of the universe and the growing mode for the density contrast 
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in order to discern among more suitable cosmological model, it is 
given by: 



0+F 

A cosmological model with better characteristics to represent 
the universe must be satisfied if £ > I. Taking into account the 
growing mode (31) and relation (33), the condition (40) implies that 
(3 < 0. Negative values for the parameter [3 violates the second law 
of thermodynamics. Then we conclude that the respective cosmo- 
logical model to the source given by ^ = 3nj3H does not give a better 
result than the usual FRW model when the age problem and the 
reach of the non-linear regime are taking together. On the another 
hand if we examine the difference of the age for the universe and 
the reach for the non linear regime in the context of the open sys- 
tem cosmology with a particle source given by ^ = 3n/3H, we note 
that this difference is smallest than the result obtained in the FRW 
framework. 
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